Exploring the breakup and transfer coupling effects in ^Be elastic scattering 
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Two cluster structures of ''Be namely, ^He+''He and *Be+n have been considered for its breakup 



to describe the available data on elastic scattering for ^Be projectile with ^^Si, 



*Zn and ""Sm 



targets. The results of these calculations suggest that the breakup coupling effects are significant 
for ''He+^He cluster model at above the barrier energies, while they are dominant at relatively 
lower energies for *Be+n model. The addition of one neutron stripping channel in the *Be+n model 
gives an overall good description of the elastic data for all the systems considered. The couplings 
generated by breakup in *Be+n model have different behaviour than the coupling effects obtained 
within the ^He+''He model, the former being more prominent at lower energies and for the heavier 
target systems. The behaviour of extracted dynamic polarisation potentials (DPP) generated due 
to breakup and one neutron transfer couplings have been investigated. 



PACS numbers: 25.70.Bc, 24.10.Eq, 21.60.Gx, 25.70.Hi 



I. INTRODUCTION 

The study of reaction dynamics for the weakly bound 
nuclei offer ample opportunities to understand the un- 
derlying cluster structure effects. The measurements of 
elastic scattering and breakup observable in the experi- 
ments involving such nuclei show that the scattering and 
breakup probabilities are sensitive to the internal clus- 
ter structure. The reactions induced by stable weakly 
bound nuclei ^'^Li, have been explained quite success- 
fully in terms of their predominant a -I- d and a + t 
cluster structures, respectively However, the clus- 

ter structure of ^Be which can be thought in terms of two 
a particles and a neutron is not clear. Within the three- 
body a -|- a -|- n picture of ^Be, no two constituents alone 
can form a bound system, a case analogous to that ob- 
served in some loosely bound unstable nuclei such as ^Hc 
and ^^Li, the so called borromean structure. The clus- 
ter structure of ^Be is not only of interest for studying 
the reaction dynamics of loosely bound nuclei but is of 
relevance to certain aspects of nuclear astrophysics. The 
nucleosynthesis via the reaction path a-t-a + n— s-^Be 
followed by ^Bc(a, n)^^C 01 is the most efficient path to 
bridge the stability gaps at A = 5 and A = 8 and knowl- 
edge of cluster structure of ^Be is vital for calculating 
the reaction rates. While a three body a -I- a -I- n cluster 
structure picture of the ^Be nucleus is more accurate one, 
the effective two-body a + ^He [a -I- (a -I- n)] d, @ or 
n -I- *Be [n -I- (a -I- a)] [HI cluster configuration can 
be used successfully in the calculations to explain the re- 
action dynamics. Experimentally, several measurements 
of ^Be breakup jl2Hl5| have been performed to quantify 
the contribution of different cluster decay components in 
its low-lying excitation spectrum. 

There are many experimental evidences which show 
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that the reaction dynamics of weakly bound nuclei are 
quite different from the tightly bound nuclei. There is 
a greater importance of the direct reactions, such as 
breakup or transfer, which may be enhanced owing to the 
low binding energies and the favourable Q values for se- 
lected transfer channels. The coupling effects due to the 
low-lying resonances and non-resonant continuum aris- 
ing due to small breakup threshold of a weakly bound 
nucleus arc expected to be quite dominant. In addition, 
the coupling effect due to transfer channels are also found 
to be significant. The Continuum Discretized Coupled 
Channel (CDCC) formalism presents an effective method 
to take into account the coupling effects of breakup on 
elastic scattering and fusion process. In this method, 
the coupled channels calculations arc performed includ- 
ing the bound states of projectile and the continuum of 
its excitation which is discretized into a finite number 
of bins. Further, the transfer mechanism can be studied 
through the Coupled Reaction Channels (CRC) formal- 
ism. The coupling effects are manifested in the behaviour 
of the equivalent Dynamic Polarization Potential (DPP), 
which may provide a qualitative idea of major couplings 
that have significant effect on the reaction dynamics [1^1 ■ 
The nature of DPP arising due to the couplings of the 
breakup and transfer processes is not clear in case of 
weakly bound nuclei. In general, the breakup couplings 
are assumed to give repulsive behaviour for the real part 
of DPP. However, the opposite effect is observed at low 
energies where the breakup couplings have the dominant 
dipole contribution. Similarly, while the transfer cou- 
plings are expected to give attractive couplings, but in 
case of positive Q-value reactions, repulsive couplings 
have been found • In addition, it has been shown 

that the couplings due to breakup also change their na- 
ture as a function of energy [l6l [l8| . 

In a recent study, we have thoroughly investigated 
two cluster configurations of ^Be namely, the 
^He-h^He (breakup threshold = 2.46 MeV) and the 
®Be-|-n (breakup threshold = 1.67 MeV) through high 
precision elastic scattering data at sub-Coulomb barrier 
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energies on ^°®Pb target. It was demonstrated that, at 
these sub-Coulomb barrier only ®Be-|-n cluster structure 
of ^Be is able to explain the data satisfactorily. Also, 
the effect of one neutron transfer was considered in these 
calculations and found to be dominant even at such sub- 
Coulomb barrier energies. It is to be noted that, in the 
previous study by Keeley et al. @, the cluster structure 
of ^Hc-|-^He for ^Be was used to explain the clastic scat- 
tering angular distribution data of the ^Be-|- ^ospj^ gyg_ 
tem, at energies around and above the Coulomb barrier. 
It was also pointed out that the coupling effect of one 
neutron stripping channel plays an important role. The 
competing coupling effects of one neutron transfer and 
the breakup channel for the ^Be projectile is crucial for 
the proper understanding of ^Be elastic scattering data 
for different target systems. Motivated by these investi- 
gations, we have attempted to validate the two models of 
9Be for few other targets (^^Si ^"^Zn H and "''Sm 
[U) in different mass regions, where elastic scattering 
data is available around Coulomb barrier energies. 



II. CALCULATIONS 

To study the effect of breakup and transfer couplings, 
we have carried out detailed coupled channels calcula- 
tions. The breakup of the ^Be in the reactions for the 
^Be-|-^*Si, ^''Zn and ^*^Sm systems has been taken into 
account by performing the CDCC calculations. In addi- 
tion to breakup (BU) couplings, the effect of one neutron 
transfer couplings have also been investigated through 
the Coupled Reaction Channel (CRC) calculations. The 
code FRESCO version FRXY.h [H is used for these cal- 
culations. The CDCC calculations are performed con- 
sidering ^Be as ^He-l-^He and *Be-|-n clusters. In the 
^He-f^Hc cluster picture of ^Be, the ground state of ^Be 
(3/2 ) is constructed by taking the relative angular mo- 
mentum L = and 2 between the core ^He (3/2~) and 
the ^He (0+) cluster. The L = 2 component is taken in 
order to account for the reorientation of the highly de- 
formed ^Be nucleus. The bound state and the resonances 
are generated by using the potential between ^Hc (3/2~) 
and ^He (0+) clusters, taken from Ref. [1]. The spectro- 
scopic amplitudes 0.81 and 0.5358, obtained from a shell 
model calculation arc used for the two components L = 
and L = 2 of the g.s. wave function, respectively Q. Sep- 
arate g.s. wave functions are obtained for the L = and 
L = 2 components by adjusting the binding potentials 
of '^Hc + ^He configuration to reproduce the experimen- 
tal binding energy of ^Bc. While the deformation in the 
g.s. is accounted by taking both the L = and L — 
2 components explicitly, the continuum states are calcu- 
lated as purely L = 0, 1, 2 states. This approximation 
which omits some terms of the full orthogonal combina- 
tions is seemingly non-trivial. However, in addition to the 
ground state, the 5/2~ inelastic state at energy 2.43 MeV 
and the 7/2~ resonance state at 6.38 MeV are generated 
using the ^He-|-"'lle cluster model. The breakup calcu- 



lations including these states along with the ^Hc+^He 
non-resonant continuum are performed. The ^He+^He 
continuum model space in momentum is limited to < k 
< 0.8 fm~^. The discretization scheme is suitably modi- 
fied to take into account the resonances. 

In the alternate cluster picture, the ^Be-|-n cluster 
structure of ^Be is assumed. The ground state wave- 
function of ^Bc is generated by coupling the valence neu- 
tron in lp3/2 state to the ^Be (0+) core configuration. 
The Woods-Saxon potential parameters (radius and dif- 
fuseness) along with a spin-orbit component for the bind- 
ing of neutron in ^Be is taken from Ref. [2^. The 1/2^ 
and 5/2'*' resonance states arc generated by using the 
radius and diffuscness parameters same as that of the 
ground state while the potential depth is varied. The 
non-resonant continuum states are generated using the 
same potential as that of the resonance states. The fi- 
nal CDCC calculations are performed by including the 
non-resonant continuum and the resonance states. The 
cluster folding (core-target and valence-target) potentials 
required in CDCC calculations for constructing ^Be -I- 
targct interaction potential arc taken from Rcfs. [2l|, [23 - 
[29I as given in Table ID In the final calculations, the depth 
of real part of optical potential for n-|-^^'*Sm is normalised 
by a factor 0.8 at ^Be incident energies 37 MeV and be- 
low, which is needed to explain the data satisfactorily. 
The same factor is also used for re-normalizing the depth 
of real part of optical potential in n-f-^^Si at the lowest 
(12 MeV) energy. 

In addition to CDCC calculations for breakup, the 
CRC calculations for the one neutron stripping channel 
are simultaneously performed to study its effect on elas- 
tic scattering for the ^Be-|-n model. The CDCC wave 
function calculated as described previously is used in the 
post-form transfer transition amplitude [30l [3l| to in- 
clude coupling of the BU states to the transfer channels. 
The optical model potentials used in the exit channels 
are same as the *Be -I- target potential parameters as 
listed in Table U in both the cases. The neutron strip- 
ping channel viz; 28Si(9Be,8Be)29Si, ^^Zn{^Be,^Bef^Zn, 
"4Sm(9Be,*Be)i'^5Sm have positive Q-values 6.81 MeV, 
6.31 MeV and 5.09 MeV respectively. The excited states 
of the residual nucleus considered in the CRC calcula- 
tions arc chosen where the well defined s pec troscopic in- 
formation is available in the literature (23l . l32ll33j as listed 
in Tabic HH 

III. RESULTS AND DISCUSSION 

A. Elastic Scattering 

The measured elastic scattering data along with the 
calculations for the three systems are shown in Figs. [TJ 
[2] and [3] at a range of energies around and above the 
Coulomb barrier in each case. The calculations with- 
out couplings (uncoupled), only BU couplings, and full 
breakup plus one neutron transfer (BU-TR) couplings 
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TABLE I: Optical model potentials used in the CDCC and CRC calculations. 



System 


Vo 


ro ao 


Wo 


ro ao 




la 


as Ws rs 




Ref. 




(MeV) (fm) (fm) (MeV) (fm) (fm) (MeV) (fm) 


(fm) (MeV) (fm) 


(fm) 




*He+2*Si 


50.7 


1.25 0.81 


20.7 


1.63 0.51 










[28j 


*Be+2*Si 


40.0 


0.89 0.87 


72.7 


0.89 0.87 










[29] 


28 c- 

n+ bi 


141.5 


1.29 0.60 










2.1 1.29 


0.60 


[27| 


''He+''''Zn 


113.6 


1.57 0.46 


15.0 


1.67 0.18 










m 


SBe+'^^Zn 


126.0 


1.10 0.60 


17.3 


1.20 0.75 










[25] 


n+s^Zn 


70.0 


1.28 0.57 






2.5 


1.28 


0.57 6.2 1.28 


0.57 


[27] 


''He+i4-*Sm 


50.5 


1.47 0.59 


18.7 


1.49 0.65 










[26j 


8Be+"*Sm 


140.0 


1.06 0.71 


112.0 


1.06 0.71 










[2^ 


n+"*Sm 


70.0 


1.30 0.58 






2.2 


1.30 


0.58 5.5 1.26 


0.58 


mi 



a 

"For ^He+targot, the same potential of '*He+target was used with 
the diffuscncss parameter increased by 0.1 fm 



TABLE IL Energy levels of residual nuclei and spectroscopic 
amplitudes (SA) used in the CRC calculations. 



^«Si 


"''Zn 


^*^Sm 


E 


J" 


SA 


E 


J" 


SA 


E 


J" 


SA 


(MeV) 






(MeV) 






(MeV) 






0.00 


1/2+ 


0.69 


0.00 


5/2^ 


0.46 


0.00 


7/2- 


0.78 


1.27 


3/2+ 


0.83 


0.06 


1/2- 


0.58 


0.89 


3/2- 


0.66 


2.03 


5/2+ 


0.22 


0.12 


3/2- 


0.41 


1.11 


13/2+ 


0.81 


3.62 


7/2^ 


0.62 


0.21 


5/2- 


0.33 


1.43 


9/2- 


0.92 


4.93 


3/2" 


0.84 


0.86 


1/2- 


0.53 


1.61 


1/2- 


0.91 


5.95 


3/2~ 


0.37 


1.04 


9/2+ 


0.79 


1.66 


5/2- 


0.64 


6.38 


1/2- 


0.78 


1.35 


5/2+ 


0.50 


1.79 


9/2- 


0.58 








1.86 


1/2+ 


0.33 


2.71 


13/2+ 


0.55 








4.40 
4.78 


1/2+ 0.52 
1/2+ 0.65 









are shown with dashed, dashed-dot-dot and solid lines re- 
spectively. The effect of the breakup and transfer on the 
elastic scattering is evident from the difference between 
the uncoupled calculations and the results obtained from 
including all BU couplings and subsequently also includ- 
ing the transfer (in the ^Be-|-n model). As can be seen 
from these figures, at the above Coulomb barrier energies, 
the ^He-|-^He model gives significant breakup coupling ef- 
fects (left column of Figs. [U [2] and [3]). However, at rela- 
tively lower energies, more dominant coupling effects are 
observed in calculations with the ®Be+n breakup model 
(Right column of Figs. [U [5] and [3]). It is observed that, 
the coupling of breakup channel leads to a rise in the elas- 
tic cross sections at the backward angles for the ^He+"'He 
model. This is usually the behaviour of breakup contin- 
uum couplings for nuclear systems with weakly bound 
projectiles. In the ^Be-|-n model however, the calcula- 
tions show a reduction of elastic cross section from the re- 
gion of Coulomb rainbow to the backward angles. Within 
this model, the coupling effects of breakup is dominant 
for a wider range of energies for all the systems. In the 



test CDCC calculations we have verified that the reso- 
nances 2.43 MeV, 5/2" state in case of '"'Hc+'^IIe model 
and 1.78 MeV, 1/2+ in case of *^Bc+n model, have the 
dominant contribution in the BU coupling. 

With the inclusion of one neutron transfer channel, 
the ^Bc+n model gives an overall good description of 
the data for these three systems over the entire energy 
range. However, the transfer couplings are found to af- 
fect the elastic cross-section not as much as the BU cou- 
plings. This observation is different from the coupling 
effects observed in ^Be-|-^°®Pb system, where large cou- 
pling effects were observed at lower energy (even upto 
10 MeV below the Coulomb barrier) due to one neutron 
transfer. In the present systems, the transfer couplings 
have a general tendency to give a rise in the backward 
angle elastic cross section which is opposite to the effect 
observed in ^Be-|-^°*Pb case This can be ascribed to 
the relatively less positive Q-value for one neutron trans- 
fer (-1-2.27 MeV) which leads to better optimum Q value 
matching in the latter case. 



B. Dynamic Polarisation Potential (DPP) 

To understand the observation from the coupling ef- 
fects in the calculations for the elastic scattering angular 
distribution in a better way, we have investigated the 
behavior of the DPP generated due to these couplings. 
DPP provides a useful way to simulate the influence of 
breakup and transfer channel coupling effects by solving 
the single channel Schrodinger equation for the elastic 
scattering with an effective potential which comprises of 
the bare potential and the DPP. The real and the imagi- 
nary part of the polarization potentials generated by the 
BU-TR couplings and only BU couplings are calculated 
using the prescription of Thompson et al. [22 1. 

The calculated DPPs due to BU couplings in the 
^He-h^Hc model for the ^Bc+^*^Sm system in the vicinity 



of the strong absorption radii (Rsa) are shown in Fig. 2] 
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FIG. 1: (Color online) The comparison of measured elastic 
scattering data for ^Be+^*Si system with the coupled chan- 
nels calculations from two models: ^He+^He (left column) 
and *Be+n (right column). The dashed, dashed-dot-dot and 
solid lines are without coupling, only BU couplings, and BU- 
TR couplings respectively. 
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(Color online) Same as Fig. [T] but for ^Be+^^Zn sys- 



It is evident from Fig. |4]that the BU couplings give rise 
to repulsive real and attractive imaginary DPPs. Sim- 
ilar behaviour is also observed for the ^Be-t-^^'Zn and 
^Be-t-^^Si systems (not shown here). Next, the calculated 
DPPs due to BU and BU-TR couplings for the ^Bc+n 
model near the respective Rsa values for the ^Be-|-^''^Sm, 
^Be+^'^Zn, and ^Be-H^^Si systems are shown in Figs. [Sj [6] 
and[7]respectively. The real part of DPPs due to only BU 
couplings are found to be attractive at larger distances 
in case of ^Be-(-^^'*Sm and ^Bc-f ^^Zn systems. However, 
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FIG. 3: 

system. 
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(Color online) Same as Fig. [T] but for ''Be+"''Sm 
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FIG. 4: (Color online) Real and imaginary parts of the 
DPPs around the strong absorption radius (Rsa=11.3 fm) for 
^Be+^''*Sm system with the ""He+^He model. The potentials 
are due to the BU couplings. 



in case of ^Be+^^Si system, the real part of DPPs are 
attractive at the lowest energy (12 MeV), whereas it is 
repulsive at higher energies. 

The attractive nature of the real part of DPP ob- 
served in the case of ^Be-|-^^''Sm and ^Be-|-^^Zn systems 
is slightly reduced due to the inclusion of transfer cou- 
plings (left part of Figs. [51 [H]). For the ^Be-l-^'^Si system, 
the transfer and BU couplings considered together (left 
part of Fig. [7|) lead to an overall repulsive real DPP. 
The repulsive transfer coupling effect is also observed 
for the ^Be-l-^'^^Sm and ^Be-|-^'*Zn systems, that leads 
to an overall repulsive real DPP at interior distances. 
The repulsive coupling effects due to transfer to positive 
Q value channels has been pointed out also by Keeley et 
al. P, [l^ . To summarize the observations on DPP, we 
encounter a peculiar case where, the real part of DPP is 
attractive due to BU couplings and it is repulsive due to 
transfer couplings contrary to the conventional wisdom. 
This observation underlines the importance of studying 
the detailed nature of coupling effects for each weakly 
bound projectile on a case to case basis, before making 
predictions about the gross features observed in such re- 
actions. For example, phenomena of the fusion suppres- 
sion and breakup threshold anomaly observed in many 
weakly bound nuclei have been often ascribed to the re- 
pulsive couplings arising due to breakup mode. However, 
the repulsive couplings due to the transfer processes and 
attractive couplings due to the breakup mode as observed 
here may call for a closer look on these effects before mak- 
ing such inferences or making any further predictions. 
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FIG. 5: (Color online) Real and imaginary parts of the 
DPPs around the strong absorption radius (Rsa=11.3 fm) for 
^Be+^^^Sm system with the **Be-f n model. The potentials are 
due to the BU-TR couplings (left column) and BU couplings 
(right column). 
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FIG. 6: (Color onhne) Same as Fig.[5]but for ''Be+''''Zn sys- 
tem around Rsa=9.4 fm with *Be-|-n model. 



To gain further insights, the energy dependence of the 
DPP arising due to the BU and BU-TR couplings at Rsa 
has been studied. The real and imaginary parts of DPPs 
at the respective Rsa due to BU-TR couplings (solid cir- 
cles) and only BU couplings (empty circles) for ^Be-f^^Si, 
^Be-|-^^Zn and ^Be-|-^*^Sm with *Be-|-n model are shown 
in Fig. [5] These calculations have been done with the 
energy independent potentials listed in Table |I] without 
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FIG. 7: (Color online) Same as Fig.[5]but for ^Be+^^'Si system 
around Rsa=8.2 fm with *Be-fn model. 



any re-normalization at any energy. From these figures, 
it is again evident that the inclusion of transfer couplings 
reduces the strength of real DPPs. The real DPP for the 
lighter system ^Bc-|-^^Si, is found to be repulsive while 
for the ^Be-|-^^^Sm system it is attractive at all energies. 
The attractive DPP in the latter case implies that a usual 
threshold anomaly could be present for the ^Be-l-^^^'^Sm 
system, i.e. at lower energies the real part of the total po- 
tential increases in strength (becomes more attractive), 
while the imaginary part decreases at energies around 
the Coulomb barrier. However, from the Fig. [5] it is ob- 
served that the imaginary part of the potential shows an 
increase in strength at energies near the Coulomb barrier 
indicating the presence of unusual threshold anomaly. A 
similar behaviour in the effective potential has been ob- 
served by the optical model analysis of measured data for 
the 9Bc4-^°^Pb,209Bi systems [SJ. The additional cou- 
plings arising due to transfer lead to an effective DPP, 
the real part of which is less attractive compared to that 
given by only BU couplings. While the real part of DPP 
given by the BU-TR couplings still remains attractive for 
the relatively heavy target system ^Be-l-^^^'^Sm, it turns 
more repulsive for the light target system ^Be-|-^^Si. The 
^Be-|-^*Zn system is an intermediate case, where the real 
DPP at Rsa is small attractive due to only BU couplings 
and turns repulsive due to inclusion of transfer couplings. 
Therefore, it seems that there is a continuous evolution 
in the behaviour of the real part of DPP from the attrac- 
tive real for the ^Be-f ^'*'*Sm system, to the case of light 
target system ^Be-|-^®Si, where this is repulsive. Here 
we would like to remark that the extent of increase of 
the attractive real potential may be somewhat reduced 
if the contributions due to continuum couplings due to 
^Hc-f-'*He model arc also significant. 
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FIG. 8: Real and imaginary parts of DPPs at the respective Rsa due to BU-TR couplings (solid circles) and only BU couplings 
(empty circles) for ^Be+'^^Si, ^Be+^^Zn and ^Be+^**Sm systems with *Be+n model. The solid lines are used to guide an eye. 



IV. SUMMARY 

To summarize, the effect of breakup and transfer cou- 
plings have been studied on elastic scattering in the 
^Be+^^Si, ^^Zn and ^'^''Sm systems. The elastic scatter- 
ing data available around Coulomb barrier have been uti- 
lized for these investigations. Two different cluster mod- 
els of ^Bc, namely, ^He-|-'*He and **Bc+n have been inves- 
tigated. CDCC (breakup) calculations for both the mod- 
els and CDCC-CRC calculations for the breakup plus 
one neutron transfer in the ®Be+n model have been per- 
formed. We obtained a good agreement of the coupled 
channels calculations with the available data. The *Be-|-n 
model gives an overall good description of the data over 
the entire energy range for all the systems. The calcula- 
tions for the ^Be-|-n model show significant BU coupling 
effects at lower energies, whereas for the ^He-|-''He model, 
large BU coupling effects are seen at energies above the 
barrier. Detailed investigations using the dynamic polar- 
ization potentials show that the attractive real DPP is 
obtained for the *Be-|-n breakup model for all the sys- 
tems at energies below the barrier while the repulsive 
real DPP is obtained for the breakup via the ^He+'*Hc 
model. In general, larger attractive real DPP due to cou- 
plings of the ®Be-|-n breakup are obtained in more heavier 
systems which remains attractive even at energies above 
the barrier for these systems. Inclusion of coupling to the 
single neutron transfer channel in the ^Be-|-n model also 



gives interesting effects in the present study. In the com- 
bined CDCC-CRC calculations, repulsive real DPP are 
obtained due to transfer processes at all energies, while 
an overall attractive real DPP is generated due to the 
dominance of the breakup in ^Be+n model at lower en- 
ergies. However, the imaginary part of the DPP shows 
an increasing trend towards the lower energy consistent 
with the phenomenon of unusual threshold anomaly. 

An extension of the present work could be the com- 
parison of present results with those obtained from the 
calculations using the three body model of ^ Be in a full 
four-body CDCC calculation. Such a calculation could 
provide a unified framework for understanding the effects 
due to both the ^He+'^He and *Be+n models apart from 
studying the possible pure three body effects. 
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